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Our work originated from harmonic analysis on the infinite symmetric group. The 
problem of spectral decomposition for certain representations of this group leads to 
a family of probability measures on an infinite-dimensional simplex, which is a kind 
of dual object for the infinite symmetric group. 

To understand the nature of these measures we interpret them as stochastic point 
processes on the punctured real line and compute their correlation functions. 

The correlation functions are given by multidimensional integrals which can be 
expressed in terms of a multivariate hypergeometric series (the Lauricella function of 
type B). 

It turns out that after a slight modification ('lifting') of the processes the correla- 
tion functions take a common in Random Matrix Theory (RMT) determinantal form 
with a certain kernel. 

The kernel is expressed through the classical Whittaker functions. It depends on 
two parameters and admits a variety of degenerations. They include the well-known 
in RMT sine and Bessel kernels as well as some other Bessel— type kernels which, to 
our best knowledge, are new. 

The explicit knowledge of the correlation functions enables us to derive a number 
of conclusions about the initial probability measures. 

We also study the structure of our kernel; this finally leads to a constructive 
description of the initial measures. 

We believe that this work provides a new promising connection between RMT 
and Representation Theory. 



In this paper we review our results from [Part I - Part V] . 

1. The Thoma simplex [Part I, §1]. The starting point of our study was the 
work [KOV] on generalized regular representations of the infinite symmetric group. 
These representations depend on two parameters (further denoted by z and z'). 
The decomposition of the representations into irreducibles is governed by certain 
probability measures P zz ' living on the infinite-dimensional simplex 

oo 

fl = { ai > a 2 > . . . > 0; ft > ft > . . . > | ^(a, + ft) < 1} 
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called the Thoma simplex [VK], [KV]. Note that O is compact in the topology of 
pointwise convergence. 

Our aim is to understand these measures. Our results show that the measures 
P zz i are close to stochastic point processes arising in scaling limits of random matrix 
ensembles. 

2. Probability measures P zz > [Part I, §§1,2]. The measures P zz i are defined 
as follows. There exists a family {s\} of continuous functions on fi, indexed by 
the Young diagrams A. Their linear span is dense in C(O). We know explicitly the 
integrals 

fzz'W = / s\(u)P zz >(du) (2.1) 
Jo, 

which determine P zz > uniquely. Essentially, this is the only information about P zz > 
that we possess. 

Now we shall describe the functions s\ and write down the formula for <p zz '. 
Let <jj = (a |/3) range over O. Set 



1, k = 1 

oo oo 

£«? + (-i)* -1 x:#, k>2 

i=l i=l 



and for any partition p = (pi, . . . , pi) 

These are continuous functions on O. The functions are related to p p 's in exactly 
the same way as the Schur functions s\ are related to the power sums p p , see [M, 
§1.7]. Specifically, 

p 

where p ranges over the partitions of |A| = Ai + A 2 + . . . ; x X is the irreducible 
character of the symmetric group of degree |A|, Xp is hs value on the conjugacy 
class indexed by p, and z~ x ■ |A|! is the cardinality of this conjugacy class. 
Conversely, 

p P = J24^- ( 2 - 2 ) 

A 

The functions ls\ are called extended Schur functions, see [KV] . 

To define <p 22 '(A) we need Frobenius notation for a Young diagram A: 

A = (pi,... ,Pd\qi,--- ,qd); 
here d is the number of diagonal boxes of A, and 

Pi = Xi- i, q i = \' i - i, 
where A' stands for the transposed diagram (see [M, §1.1]). 
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Set n = |A|, t = zz'. Then 

V - W= W,U (23) 

Ui<j[(Pi ~Pj)(Qi ~Qj)] 
YlijiPi + Qj + l) 

where (a) m = a(a + 1) • • • (a + m — 1), (a)o = 1. 

We shall always assume that z and z' satisfy one of the following conditions 

(*) z' = z, z e C\Z; 

, , ( 2 - 4 ) 

(**) z and 2 are real and for a certain m G Z, m < z, z < m + 1. 

Under this assumption there exists a probability measure P zz > such that (2.1) holds. 

Note that the measure P zz > can be obtained as the limit, as n — > oo, of certain 
statistics on partitions of n, see [KOO, proof of Theorem B], [Part II, Introduction]. 

3. Stochastic point processes [Part I, §4]. We shall interpret the measures 
P zz i as stochastic point processes on J = [-1, 1] \ {0}. With each point wGflwe 
associate a point configuration in /, 

u = (a\ 13) i-> (ai, a 2 , • • • ; -(h, -fo, • • • ) 

where we omit possible zeros in a and (3. Note that any such configuration has no 
accumulation points in /, the points condensate near the origin which is not in I. 

Thus, the measure P zz / becomes a probability measure on the space of point 
configurations in /, i.e., a stochastic point process on /. We shall denote this 
process by V zz > . 

The nth correlation function p n (xi, . . . ,x n ) of a point process is the density of 
the probability to find a point in each of the infinitesimal intervals [xi,Xi + dxi\ 
(see [DVJ], [Part I, §4] for details). 

Our aim is to compute the correlation functions pn Z of our processes V zz '. 

4. Moment problems [Part I, §3], [Part II, Section 1.1]. Our strategy to 

( zz' ) 

solve the problem is to find the probability measures a\ on [— 1, l] n , n = 1, 2, . . . , 
characterized by their moments 

f Xl •■•Xn<7 ( n Z ' ) (dx) = J + ,/„+!) (u)P ez > {du>) . 
[-1,1]" Q 

Note that the RHS is explicitly known because of (2.1), (2.2): 

/ x^...x l -a^'\dx)= J2 4 1+ i,...,/ n+ i)P«'(A). (4-1) 
[-1,1]" A - |A|=^)(ii+i) 

i 

It turns out that outside the diagonals 

p^'Hxx, ... ,x n ) = r {density of a^ z '\dx u . . . ,dx n )} (4.2) 



\X\ ■ ■ ■ X 



n | 
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( zz') 

It is worth noting that the measure a n ' for n > 2 always has singular components 
living on the diagonals while (as is proved in [Part I, Proposition 4.2], [Part II, 

Theorem 2.5.1]) p£ z ^ has no such components 1 . 

Thus, we obtain the correlation functions from a more sophisticated object. 

(zz') 

However, we can not formulate a moment problem for p n directly because of the 
absolute value sign in (4.2). 

5. Integral represenations for the correlation functions [Part I, §§5,6], 
[Part II, Chapter 2]. The moment problem (4.1) can be completely solved. In 

particular, for n = 1 we obtained a two dimensional integral represenation of a[ zz ^ 

(see [Part I, Theorem 5.8]). A careful examination shows that a[ zz ^ has no atom 
at zero. This fact has an important corollary. 

Theorem I ([Part I, Theorem 6.1]). With probability 1, 

oo 

+ A) = L 

i=l 

For n > 1 the solution of the moment problem (4.1) requires a lot of work ([Part 
II, Chapter 1]). The difficulties come from the fact that the RHS of (4.1) involves 
the symmetric group characters for which there is no simple expression. We use 
Murnaghan-Nakayama rule ([M, §1.7, Ex. 5]) to handle the characters. In [Part 
II, Theorem 1.2.1'] we obtain a more economic form of this rule which enable us to 
solve the moment problem. 

(zz') 

The final expression for p n is a linear combination of multidimensional inte- 
grals of various orders up to 3n. The situation simplifies when all the variables 
xi, . . . ,x n are of the same sign, say, positive. 

Theorem II ([Part II, Theorem 2.2.1]). Let xi, . . . , x n > 0. Then 

p n zz '\x 1 ,...,x n ) = t n T(t) 

J fj. r (_ z + i) T (z' + 1) T(-z' + 1) T(z + 1) 

ai,bi>0 ' (5.1) 



X 



^x i (a i +b i + l)<l 

x det | 1 ^ (i-E.»^ + VD)'-- dadk 



ai + bj + 1 / T(t — n) 

The RHS of (5.1) is well defined for z and z' such that 

-1 < Uz, $lz' < 1; t = zz' > n. 

For other values of z,z' we use analytic continuation. 

— t 1 \ 

We tacitely assume that Y2i x i < 1- Actually, the correlation measure pn (x)dx 

always lives on the set Yli \ x i\ — 1> additional considerations show that there are 

no singular components on the faces J2i \ x i\ = 1) see the beginning of the proof of 

Theorem 3.3.1 in [Part II]. 



1 This just means that nonzero coordinates in a and (3 are pairwise distinct with probability 1. 
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The RHS of (5.1) can be expressed via the multivariate Lauricella hypergeometric 

function of type B, see [Part II, Section 2.4]. In particular, p\ can be expressed 
in terms of the Appell's two-dimensional hypergeometric function F3, see [Part II, 

(zz') 

Corollary 2.4.2]. Another expression of p\ through the Lauricella function in 
three variables is given in [Part I, Theorem 5.12]. 

6. Lifting [Part II, Chapter 3]. A surprising fact is that the correlation func- 
tions are greatly simplified after a 'lifting' of our processes V zz i to a slightly bigger 
space. 

According to Theorem I, the measure P zz * ic concentrated on the face 



on Qq. In other words, we tensor P zz > by the gamma-distribution with parameter 
t. 

We associate to a point u = (u>, s) G O a point configuration in R* = R \ {0} 
as follows 



Thus, we get a probability measure on the space of point configurations in R*, i.e., 
a point process. We denote it by V zz ' ■ 

The process V zz i is obtained from V zz i by multiplying the random configuration 
in / by the scalar random factor s with gamma-distribution. 

The lifting is invertible via the map 



where s = J2i( a i + A)- 

Note that s is finite almost surely with respect to P zz i . 

Let p^ z \xi, . . . , x n ) be the nth correlation function of V zz ' ■ 

The 'lifted' correlation functions are related to the initial ones by the following 

simple transformation 



n = {u = (a\/3)\ + Pi) = 1} c a 



Let 



n = O x R+. 



We pass from the measure P zz > on Oq to the measure 




(H P),s) i->(ai, a 2 , -/3 2 , • • • ) 

= (sai, sa 2 , . . . ; -s/?i, -s/? 2 , ■ • • )• 




oc 




(6.1) 
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see [Part II, Proposition 3.1.1]. This implies that the moments of the correlation 
measures ph* \x)dx and ~pfi z \x)dx are related in a very simple way 

J x^ ■ ■ ■ x™ p^ ^ (xi , . . . , x n )dxi . . . dx n 

= (t)i 1+ ...+i n J x 1 ? ■■■x l ^p^ z '\x 1 ,... ,x n )dx 1 ...dx n 

(the moments are finite if l±, . . . , l n > 1). 

Note also that the transformation (6.1) is invertible, and there exists a complex 
inversion formula similar to that for the Laplace transform. 

7. Determinantal formula. Matrix Whittaker kernel [Part IV, Sections 

~ (zz r ) 

1,2]. Now we shall show that the correlation functions p\ ' have determinantal 
form with a kernel expressed through the Whittaker function W KjAt (x), x > 0. 
This function can be characterized as the only solution of the Whittaker equation 

'-(H+^)"» 

such that y ~ x K e~% as x — > +oo (see [El, Chapter 6]). Here n and p are complex 
parameters. Note that 



We shall employ the Whittaker function for real k and real or pure imaginary p; 
3n VI 
Set 



then W Kjfl is real 



z + z' - - : 



a= — - — , p = — - — , a = ^sm(nz) sm(irz') > 0. (7.1) 

It is often convenient to consider (a, p) as new parameters of our processes, instead 
of (z, z'). In terms of (a, p) the basic assumptions (2.4) take the following form 

• a is always real 

• either p is pure imaginary 

• or p is real, \p\ < -, and there exists m G Z such that 

m + \p\ < a < m + 1 — \p\ 

• a is not an integer when p = 0. 

Introduce the following functions on R + 

A+(x) = y/xW a+ i jlt (x) A-(x) = y/xW_ a+ i tli (x) 

B+(x) = V^W a _iJx) B_(x) = ^W_ a _ijx) 

Theorem III ([Part IV, Theorem 2.7]). The correlation functions of the lifted 
process V zz ' have the form 

p^ z '\u u ... ,u n ) = &et[K(u i ,u j )] n ij=1 , 

n = 1,2, . . . ; it!, . . . ,u n e R*, 
6 



where the kernel K(u, v) is conveniently written in the block form 



K(u, v) = < 



' K++(u,v), 
K + _(u, -v), 
K_+(-u,v), 
I K—(-u,-v), 



u, v > 0; 
u > , v < 0; 
u < , v > 0; 
u, v < 0; 



with 



1 A + (x) J B + (y)- J B + (x)A + (y) 



a A + (x)A_(y) + tS + (a;)5_(y) 
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x + y 



K _ ( x y ) - ° A +(V) A -( X ) + tB+(y)B.(x) 
+ n x + y 



(7.3) 



K—(x,y) 
(recall that t = zz' ). 



1 



T{-z)T{-z>) 



A_ {x)B_{y)-B_{x)A_{y) 
x-y 



This is one of our main results. 

The matrix representation of the kernel 



K = 



K++ K + _ 



(7.4) 



correponds to the splitting M* = R + U R_ of the phase space and subsequent 
identification of M_ with the second copy of 1R+. We call (7.4) the matrix Whittaker 
kernel. 

All the blocks of (7.4) are real valued kernels on R + . They possess the following 
symmetry 

K++(x,y) = K++(y,x), K__(x,y) = K__(y,x), 
K+-(x,y) = -K-+(y,x). 

Note the minus sign in the last relation. It means that the kernel (7.4) is formally 
J-symmetric where J is the operator id © (-id) in L 2 (IR_|_, dx) © L 2 (R + , dx). 

Another symmetry property: the change of parameters (z,z') — > (—z,—z') is 
equivalent to the transform of the kernel 



K++ K + . 



-K-+ 



Determinantal form for the correlation functions (like (7.2)) appears in different 
problems of random matrix theory and mathematical physics, see, e.g., [Dy], [Mel], 
[Mai], [Ma2], [TW1-3], [L], [KBI]. In most situations the kernel K is symmetric 
or Hermitian (see, however, [B]). Appearance of J-symmetric kernels seems to be 
new. 
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8. The L— kernel [Part V, §2]. Consider the operator K in the Hilbert space 

L 2 (R*,du) ~ L 2 (R + , dx) © L 2 (R + , dx) (8.1) 
given by the kernel (7.4). 

Theorem IV ([Part V, Theorem 2.4]). Assume \a\ < \. Then 

K =ih <8 ' 2 » 

where L is bounded and given by the kernel 







(f) 



a exp (_£+£) 



(8.3) 



L(x,y)= ( x+y) * w x+y 

7T \x) x+y 

Recall that a and a were defined in (7.1). 

It is worth noting that (8.3) involves no special functions. Note also that \x = 
enters only in the scalar factor a. 

The formulas (7.2), (7.3), (8.2), (8.3) give a precise description of the process 
V zz ' and, thus, of the initial process V zz <. We consider these formulas as our main 
result. 

9. Spectral analysis [Part V, §3]. In this section we shall diagonalize the 
operators K and L. Denote by A the operator in L 2 (R + , dx) with the kernel 



^, B) = £(£)-^(z20 (9 .d 

7T \yj x + y 

It is bounded provided that \a\ < \. Let A' denote the operator in L 2 (M + , dx) with 
the transposed kernel 

A'(x,y) = A(y,x). 

Then (8.1) implies that the blocks of the matrix Whittaker kernel are expressed 
through A and A' as follows 



K++ = AA'(1 + AA')- 1 K__ = A'A(1 + A' A)- 1 
K + _ = (1 + AA')~ X A K_ + = -(1 + A A')' 1 A' [ ' 



Consider the following Sturm-Liouville differential operator on (R + , dx) depending 
on a as a parameter: 

B -=s^ + (-!)'• (9 - 3) 

We have 

V a A = AV_ a 

which implies that T> a (formally) commutes with AA' and K ++ while T>_ a com- 
mutes with A' A and K 



Consider the following functions on 



We have 



fa,m{x) = -Wa iim (x), 111 > 0. 



P«/a, m = ( a 2 + - + m 2 ) f a> 



According to [W], the functions f a ,m with a fixed form a continual basis in L 2 (M + , <ix)| 
diagonalizing T> a . Moreover, the Plancherel formula looks as follows. For good 
enough functions f(x) and g(x) 



oo 



(i — im ) 1 ( -—a — im\ dm 



where (•, •) stands for the inner product in L 2 (R + ,dx). 
We have 

Af a ,m = - T ( - - a + im) T ( - - a - im ) f a , m , 
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A'f a ,m = ^ T Q + a + irnj V Q + a - irnj /_ a>m . 

Returning to the decomposition (8.1), we take {/ a ,m}m>o as a basis in the first 
summand L 2 (R + ,dx) and {/_ ajm } as a basis in the second one. 

Then we get a basis in the whole space L 2 (R*,du) diagonalizing both 



A 
-A' 



and K 



AA'il + AA')- 1 (l + AA'^A 
-(l + A'A)- 1 ^ A'Ail + A'A)- 1 



In particular, we get the diagonalization of K ++ : 

cos(2n/j,) — cos(27ra) 
cos(27r/i) + cos(27um) 

Note that for a = the integral transform inverse to / 1— > {(/, f a ,m)}m>o is the 
Kontorovich-Lebedev transform, see [E2]. 

The above results show that the operators K and L with a fixed and fx varying 
form a commutative family. 

10. Applications [Part III, Sections 2,5]. Now we shall give applications of 
the main results. These applications concern our initial object — the probability 
measures P zz i on Cl. 

Theorem V ([Part III, Theorem 5.1]). Consider a±, a.2, ■ . . ; Pi, 02, • • • o,s ran- 
dom variables on the probability space (Cl, P zz ')- Then with probability 1 there exist 
the limits 

k l ™ a * =£sa =ex H — 2^)- (10J) 

This result is a kind of the strong law of large numbers. Roughly speaking, 
it means that both ctj's and (3^s decay as geometric progressions with the same 
exponent. Similar situation occurs for the Poisson-Dirichlet process, see [VS]. 

9 



The same limit relation as (10.1) holds for cn^'s and /Vs. 

The proof of (10.1) is based on the examination of p\ zz ^ and p^ zz \ an d it 
incorporates Kingman's argument from [Ki, Section 4.2]. 

Observe that the exponent of the decay does not change under the shifts (z, z') — > 
(z + N, z' + N), N e Z. Moreover, the whole process describing the tails of the 
sequences {cii} and {{3i} turns out to be invariant with respect to these shifts. This 
periodicity is quite unexpected: initial formulas (2.3) have no apparent periodicity. 

Another application relies on the examination of p\ zz ^ alone. 
Theorem VI ([Part III, Proposition 2.2]). 



Let 



i=l 



r(-) 



Then 

E \a\ = sm(irz) sin^^') 



z + z'-l ip(-z') - i/>(-z) 



+ 



2nsm(7i(z — z')) zz' n sm(ir(z — z')) 

and E \(3\ is given by the same formula with (z, z') replaced by (—z, —z'). 

Here E is the symbol of expectation. 

11. Tail process [Part III, Sections 3,4], [Part V, §4]. Here we study the 
asymptotic behavior of the process V zz ' near the origin. The starting observation 
is that 

~{zz') I \ '' 

Pi w 



U 







\u\ 



where 



c = 



2/ja ; 



(11.1) 



(11.2) 



7T sin(27r/i) 

This asymptotic relation agrees with the rate of decay of cc^'s and /3j's, see (10.1). 

As before, we identify the phase space R* with the disjoint union of two copies 
of R_|_ and then in each copy of we make the following change of variable 



x 



£ = —clnx. 



(11.3) 



Theorem VII ([Part V, Theorem 4.1]). Let K'(£,rj) denote the matrix Whit- 
taker kernel on R obtained from the matrix Whittaker kernel (7.4) by the change of 
variable (11.3). Then there exists the limit 

lim K'(£ + M,ri + M) = £(£,77). 
Here JC(£,r]) is a translation invariant matrix kernel with /C(£,£) = 1, 



/C_ + (^) K—(Z,v) 



F(Z - V) G(£ - 17) 
-Gfa-fl F(t-rj) 



where 



(11.4) 



1 (sin(7T/x) cos(7ra)) cosh(A^) + (sin(7ra) cos(7r/z)) sinh(AC) 
^' = 2^ cosh(5C) 
and the constants A, B are as follows 



TV sin(27T£t) 
4/ja 2 



B= / 2 ^ >0, 



7rsin(2^) 
^ " 2 ^ B ~ 20-2 ■ 

We call the point process onlUR with the correlation functions given by the 
determinantal formula with the matrix kernel /C(£, rj) the tail process for V zz >. By 
the construction, it describes the behavior of 5^'s and with large % after the 
appropriate rescaling. 

In particular, the kernel 77) on R describes the tail of {5^} alone, and 

/C (£,T]) does the same for {$}. Since /C ++ = /C , the tail properties of {5j} 

and are identical. 

The same kernel /C++(£, ry) appears in the scaling limit of the correlation func- 
tions for the unlifted process V zz i restricted to (0, 1]. However, this requires more 
sophisticated considerations, see [Part II, Sections 4.2, 4.3]. 

Recall that \x is either real or pure imaginary. According to this, the constant A 
is also either real or pure imaginary. In the latter case the hyperbolic sine in the 
numerator of (11.4) turns into the ordinary sine. 

The next result is parallel to Theorem IV. 

Theorem VIII ([Part V, Proposition 4.2]). Let \a\ < \. Then 

r 

K 



1 + C 

where £ is a bounded operator in L 2 (IR, d£) © L 2 (R, d£) with the kernel 



where B is as above. 



r) a_ exp(-2aB(g-Tj)) 

2tt cosh(B(£- v )) 



a_ exp(-2aB(T ? -Q) 
2-k cosh(B(r7-£)) 



(11.5) 



Finally, note that all formulas of this section are invariant under the shifts 
(z, z') -> (z + N, z' + N), N e N (or, equivalently, (a, p) — > (a + N, /x)), cf. Section 
10. 

12. Formalism of fermion point processes [Part III, Section 1], [Part V, 
§1]. Here we shall discuss elementary general properties of the point processes with 
determinantal correlation functions [DVJ], [Mai], [Ma2]. 

Let X be a phase space with reference measure dx, K(x, y) be a kernel on X, 
and V be a point process on X with the correlation functions 

p n (x 1 , . . . , x n ) = det [Kfa, Xj)]" j=1 . 

Let 2) be a subset of X such that 



I pi(x)dx < oo. 
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Then the number of points in 2) is finite with probability 1. Moreover, the prob- 
ability to find exactly n points located in the infinitesimal volumes dxi, . . . ,dx n 
around points x\, . . . ,x n equals 

det [Lsg (xi, Xj)]™ =1 

7T n (xi, . . . ,X n )dXl . ..dx n = TTm~? \ 3 dx ± " --dXn. (12-1) 

det(l + L^J 

Here L^(x, y) is the kernel of the operator in Z/ 2 (2), dx) such that 

K _ L 9 
K % - 



1 + L 



2) 

where is the operator in L 2 (2), dx) with the kernel K(x, y) restricted to 2). 

Thus, in case of finite point configurations the operator L = j^j^ has a clear 
probabilistic meaning — it gives the distribution functions iv n . 

It is tempting to apply (12.1) in our case, when the operator L has especially 
simple form. Unfortunately, we can not do this for the whole space because the point 
configurations are infinite. Of course, we can restrict ourselves to an appropriate 
region 2), but then the simple form of our L(x, y) will be lost. 

Again, in the case of finite configurations, if the space is a disjoint union of two 
pieces and L is written in block form according to this splitting, then vanishing of 
the diagonal blocks of L (as in (8.3) and (11.5)) exactly means that the random con- 
figurations has equally many points in each of the pieces (see [Part V, Proposition 
1.7]). We do not know how to interpret such vanishing in our situation. 

13. Distribution of a.\ [Part III, Section 2]. Here we shall consider the lifted 
process V zz > restricted to M + Cl*; it is governed by the '++' block of the matrix 
Whittaker kernel, see (7.3). We call K ++ (x,y) the Whittaker kernel. 

Note that for any r > the intersection of the random configuration with 2J(r) = 
[r, +oo) is finite. By (12.1) 

Prob{5! <r}= 1 r = det(l - K^ t) ). (13.1) 

det(l + LfQ( T )) 

Here, following the notation of Section 12, the kernel of -K^t) is obtained by re- 
stricting K ++ (x, y) on [r, +oo). As was pointed out by Tracy [T2], a modification of 
the argument from [TW3, Section V.B] allows to express the Fredholm determinant 
(13.1) through the Painleve transcendent V. 

14. Comparison with random matrices. Degenerations of the Whittaker 
kernel [Part III, Sections 1,6], [Part V, §5]. There are a lot of similarities 
between the processes V zz > and point processes arising from random matrices. Ran- 
dom matrix theory leads to various kernels: the sine kernel, the Bessel kernel, the 
Laguerre kernel, etc. (see [Mel], [TW1-3]) All of them have the form 

<p(x)ip(y) - ip(x)<p(y) , 14 
x - y 

and so does the Whittaker kernel K ++ (x,y), see (7.3). Here (/?(•) and ip(-) are 
solutions of certain linear second-order differential equations. The Fredholm deter- 
minants of many kernels of this form are expressed through Painleve transcendents, 
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see [TW3, Sections III and V]. Various kernels (14.1) restricted to suitable intervals 
commute with Sturm-Liouville operators (see [G], [Mel, §5.3], [TW1], [TW2]). The 
same is true for the Whittaker kernel restricted to [r, +00), see [Part III, Section 
6]. 

The Whittaker kernel degenerates to the Laguerre kernel of order iV and param- 
eter a > —1 if we formally substitute z' = N, z = N + a, see [Part III, Remark 2.4]. 
The Bessel kernel is the scaling limit of the Laguerre kernel as N — > +00. Similar 
scaling procedure applied to the Whittaker kernel leads to a two-parametric family 
of Bessel-type kernels, see [Part V, Theorem 5.1]. 

Likewise, the stationary kernels from Theorem VII generalize the sine kernel: 
if yLt — > zoo, (11.4) tends to sin 7 [^- > irrespective to the behavior of a. The kernels 
/C_i__i_(£, 77) have already appeared in [BCM, MCIN] in connection with so-called 
g-Hermite ensemble. 

As for the matrix Whittaker kernel, it has a similarity with matrix kernels arising 
from two-matrix ensembles ([EM], [Me2], [MS]), see [Part IV, Section 3]. 

Acknowledgements. We would like to thank Craig A. Tracy for his letter [T2] 
and for drawing our attention to the works [MTW] and [Tl]. 
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